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Abstract 

We consider the possibility of getting accelerated expansion and w — — 1 crossing in 
the context of a braneworld cosmological setup, endowed with a bulk energy-momentum 
tensor. For a given ansatz of the bulk content, we demonstrate that the bulk pressures 
dominate the dynamics at late times and can lead to accelerated expansion. We also 
analyze the constraints under which we can get a realistic profile for the effective equation 
of state and conclude that matter in the bulk has the effect of dark energy on the brane. 
Furthermore, we show that it is possible to simulate the behavior of a Chaplygin gas 
using non-exotic bulk matter. 



Recent experimental results suggesting an accelerating expansion of the universepQ, 
[2], have become the subject of active research. During the last few years, a number of 
theories were proposed to answer the question of the origin of such an acceleration. The 
most common explanation is to assume the existence of some, for the moment unknown, 
"dark energy" component, which drives the expansion. The physics of this dark energy 
is still speculative, as no direct experiment can provide us with actual information about 
its nature. We can only indirectly probe its properties by studying the dynamics of 
cosmic expansion. There are several theories which attempt to describe dark energy. 
Among them are the phantom [3j and quintessence [I] theories, quiessence (cosmological 
constant) and other forms of exotic matter [5], like the Chaplygin gas [6] (see also [7] 
for a review and a more complete list of references). In these theories, dark energy is 
attributed to the presence of new fields with non-trivially modified equations of state. 

One can follow a different route and try instead modified theories of gravity. These 
are theories which extend or modify GR in such a way, so that dark energy is not 
some new -maybe exotic- matter constituent, but merely a manifestation of the modified 
dynamics of geometry. In this category we find scalar-tensor theories of gravity, the so- 
called modified gravity models and theories with extra dimensions [8], [9]. Braneworld 
models, where our four- dimensional world is realized as a brane embedded in a higher 
dimensional space (bulk), have received considerable attention mainly due to the hope of 
tackling hierarchy problems and explaining the large disparity between the Planck and 
the electroweak scale pU]. Brane models like these are also directly connected to String 
Theory. An interesting consequence of brane models is that they allow for the presence 
of five-dimensional matter which propagates in the bulk space and may interact with the 
matter content of the braneworld [12J. Such an interaction can alter the profile of the 
cosmic expansion and lead to a behavior that would resemble that of dark energy. 

It has already been shown [13J that a configuration of bulk matter resembling a per- 
fect fluid moving along the fifth dimension is capable of producing accelerated expansion 
and give an effective equation of state which can cross the w — — 1 line. However, for the 
later to occur, the existence of a large and negative dark radiation term was needed. The 
presence of such a term may cause problems in early times, as it could lead to a nega- 
tive effective energy density (although the validity of the modified Friedmann equation 
for large redshifts is in question). Our goal here is to replicate these effects by ignoring 
the dark radiation term and considering only bulk matter as the late time dominant 
term. We are also interested in determining not only under which conditions we can 
have w e ff < —1, but whether we can obtain a temporal profile of the crossing which 
is at least qualitatively in accordance with observations for < z < 1. As we will see, 
these requirements can be met, provided we have a specific ansatz for bulk matter and 
the parameters are properly set. We will also demonstrate the equivalence of this bulk 
matter to a Chaplygin gas. 

We will start by considering a braneworld cosmological model, with a 5D bulk space 
and a positive tension brane embedded in it. The action of the model is 

S = J d 5 xV^G ( 2M 3 R - A + 4 m) ) + / d 4 xV=g~ (-a + £[ m) ) . (1) 

The 5D metric is Gmn with signature (— , +, +, +, +), while is the induced metric 
on the brane and M denotes the five-dimensional Planck mass. We shall assume that the 
metric can be written as 

ds 2 = -n 2 (y,t)dt 2 + a 2 (y,t)-f ij dx i dx j + b 2 (y,t)dy 2 , (2) 
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where y denotes the fifth dimension (the brane is situated at y = 0) and 7^ is the metric 
of a maximally symmetric three-space. 

Varying the action with respect to the metric, we obtain Einstein's equations 

Qmn = Rmn — -GmnR = ^j^ 3 Tmn ■ (3) 

The energy-momentum tensor Tmn resulting from the above Action is of the form 

Tmn — ^mtv + T MN — G M n^ — 9^ & 5(y) d^N > (4) 

where T^n results from and T$ N results from . In particular, we will assume 
an energy-momentum tensor for the bulk content of the form 



T {B)M N = ( Pb^ ] , TS = ( P B a 2 1%3 _0 \ . (5) 

-n 2 P 5 P B b 2 




PbU 2 —n 2 P 5 



The quantities which are of interest here are the pressures Pb and P5, as these two enter 
the cosmological equations of motion. P5 is the term responsible for energy exchange 
between the brane and the bulk. The corresponding energy-momentum tensor for the 
brane matter is 

T(b) M = 6 M diag( _ P:P7P ^ 0) _ (6) 

Substituting the above ansatze for the metric (j2j) and for the energy-momentum 
tensor ©, © into the equations of motion ([3D, we can obtain the set of cosmological 
equations for the components of the metric |1 lj [12J [13] . These equations simplify if we 
choose Gauss normal coordinates, such that b(t,y) = 1. We can also make use of the 
freedom to take n(0,t) = no(t) = 1. After these simplifications, we get for the evolution 
equations^ on the brane (ao(t) = a(0, £)) 

P + 3(p + p)^ = -2P 5 , (7) 
a 

The first expresses energy conservation on the brane. We see that positive P5 means 
energy outflow from the brane, while for P 5 negative we get the opposite effect, with 
energy flowing from the bulk into the brane. The second equation is the Friedmann 
equation for the brane world. The "55" component of the bulk pressure appears on the 
right hand side of this equation (jSJ) affecting the cosmological evolution. 

We assume that matter on the brane has an equation of state p = wp. Introducing 
the parameters (3 = (24M 3 ) -2 and 7 = a (3 and omitting the "o" subscript from the scale 
factor, we can rewrite the cosmological evolution equation (|SJ) as 



-1, 0, 1 is the maximally symmetric internal space curvature parameter. 
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We also assume the Randall- Sundrum fine-tuning to hold, so that the effective cosmo- 
logical constant on the brane 

X ^ A+ 2^P' < 10) 
is equal to zero, i.e. A = 0. The cosmological evolution equation (EJ) can then be expressed 
as a set of two equations using a dark energy field x 

* + #-2i^) = 4 4 + i) P5+ i^- (12) 

where the defining equation for \ is 

- = - x -(3w + 2)(3p 2 - (3w + l) 7 p • (13) 
a 

The first of these equations (TTTT) is analogous to the Friedmann equation of standard 
cosmology. The dark energy field accounts for contributions to the Friedmann equation 
which do not come from ordinary brane matter. In the absence of bulk content, it is 
equal to a dark radiation term, a consequence of the presence of the extra dimension. 
When Pb or P 5 are non-zero, it also receives contributions from these terms. Thus the 
origin of \ can be both geometric and matter-related. 

We are interested in scenarios where the energy density of the brane is much lower 
than its tension, that is p « a. We can then disregard terms quadratic in p and, thus, 
we get the simplified cosmological equations 

d\ 2 k P B , 

a) = 2 ^ + *-iW' < 14 > 



L a V A 24M 3 ") ^ 3 ' 12M 3 
and, of course, the continuity equation 



p + 3p-(l + w) = -2P 5 . (16) 
a 

The functions Pb, P5 are functions of time corresponding to the values of Ps(y,t) 
and P5 (?/,£) on the brane. The energy- momentum conservation Vm^tv = cannot fully 
determine Pb and P5 and a particular model of the bulk matter is required [13]. Here 
we are going to consider a general ansatz for the bulk pressures (see also [T5] ) 

P B = Da" , P 5 = P f^j a". (17) 

where P and D are constant parameters. For this choice of pressures, the equation for \ 
can be solved exactly to give 



A a 4 (z/ + 4) Gu + 4) 
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where we have defined 5 = D/2AM 3 . Similarly, we can proceed with the integration of 
(fTtjj) and get 

C 2F 
P = ^3(TT^ - [3(1 + w) + fJ ] a ' 1 ' (19) 

Substituting (|18p and ffT§]) into the Friedmann equation (|14|) . we can write it in a con- 
ventional form as 

a\ k 87r 



a) + a*=T G » P °"- (20) 

where Gjy = S^/An = 3cr/47r(24M 3 ) 2 is the 4D Newton's constant and the effective 
energy density p e ff stands for 

C C/27 25 v 2(3w-l)F 

Mf o 3 ^) a 4 7(^ + 4) (/i + 4)[3(l + W ) + /i] 1 j 

At late times, the term with the highest power in a will dominate and we will get a scale 
2 

factor a(t) ~ t . Acceleration occurs whenever max(/x, z/) > —2. 

Using the already mentioned effective energy density, we can derive an expression 
for the effective equation of state parameter using the prescription [20j 

- -1 - l dl < 5H ^ (22) 
^ " 1 3 din a ' 

where 5H 2 = H 2 / Hq — Q m a~ 3 accounts for all terms in the Friedmann equation not 
related to the brane matter, for which we have taken w — 0. Using as our variable the 
redshift parameter z = — 1, w e ff has the general form 



1 • 
3 \-A- B(z + 1) 



where the parameters A and i? depend on the bulk pressures 

25 48DM 3 _ 2F 

= 7(^ + 4) ~ (z/ + 4)a' = (/i + 3)(/i + 4) ' ( } 

We have assumed that the dark radiation parameter C is small and that, at late times, the 
dark radiation term ~ is negligible with respect to the contribution of the pressures 
Pb and P5. Thus, we have ignored this term in our expression for to e //. 

Equation (I23p allows us to examine under which conditions we can have an effective 
equation of state which can exhibit w e ff = — 1 crossing, as is suggested by recent obser- 
vational data. In particular, we want to investigate what are the requirements that our 
model's parameters have to comply with, in order to get a qualitatively correct behavior 
for the evolution of w e ff at late times. Current observational findings indicate that at 
the moment (z = 0) we are at a stage of accelerated expansion with \q\ of order unity 
and effective equation of state w e // = —1.21. The w e ff = —1 crossing seems to have 
occurred at redshifts z = 0.2. Therefore u> e // seems to increase with z, starting with 
a value of about —1.21 at z = and crossing the w e ff = —1 line about z = 0.2. Our 
goal is to obtain constraints on the values of our parameters A, B, /1 and v and see for 
which combinations we can have a w e ff with the desired temporal profile, while having 
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accelerated expansion. The corresponding late time deceleration parameter (where we 
have neglected the contribution from brane matter as well as the dark radiation term) is 

1 h {v + 2)A + {ii + 2)B{z + l) v -» 
q H*a 2A + 2B {z + 1) 1 " 11 ' 1 ' 

By inspecting equation (I2"3"j) . we first notice that the denominator is the effective 
energy density at late times and thus it has to be positive at the interval of z we are 
interested in (in this discussion we focus on < z < 1). At z = this condition reduces 
to A + B < 0. It is, thus, obvious that A and B cannot be both positive. Taking the 
derivative of (|23l) . we obtain 

dwen = ^ ("-rf^ir 1 . (26) 

dz {A + B{z + iy-») 



All terms in the right hand side of (1261) are positive in the interval < z < 1, except 
from the product AB. In order to have an increasing w e ff, A and B must have the same 
sign, i.e. they must both be negative (A, B < 0). Next, we want w e jf to assume a value 
less than — 1 at z = 0. In order to achieve this we must have vA + fiB < 0. There are 
four possible combinations: 

1) /i, v > 0: The constraint is satisfied. 

2) fj, > 0, v < 0: vA + fiB < when J > -f . 

3) /i < 0, i/ > 0: vA + /US < when £ < - j. 

4) /i, i/ < 0: The constraint cannot be satisfied for A and 5 negative. 

So, the first three cases seem to be allowed. However, a further constraint comes from 
the fact that we want w e jf = — 1 for z = 0.2. This reduces to the relation 

^ = -§(1-2)^. (27) 

In the case (1), with both \i and v positive, the left hand side of (1271) is positive, while 
the right hand side is negative and thus it is ruled out. Cases (2) and (3) can satisfy this 
equation without providing any additional constraints for the powers \i and v. Thus we 
conclude that the only choices of parameters which are in accordance with the general 
profile for the evolution of w e // for < z < 1 are: 

1) A, B < with fi > 0, v < and - > -f . 

2) A, B < with n < 0, v > and J < -f . 

Since one of the powers has to be positive, the condition max(fi, v) > —2 is always 
satisfied and we also get the desired accelerated expansion for late times. We also see that 
the ratio of bulk pressure parameters A and B determines the ratio of the two powers. 
A large difference between the initial pressures D and F would result in an accordingly 
large disparity between fi and v. This is demonstrated in Figure 1, where we give a 
graph of the deceleration parameter (q) and w e ff, for a set of values that satisfy the 
above constraints. For this choice, we see that there is an accelerated expansion at late 
times and w e ff has a realistic temporal profile. 

The requirement that A and B must be negative has a direct impact on the sign of 
parameters D and F and the corresponding behavior of the underlying bulk content. In 
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Figure 1: Graph of the deceleration parameter q and w e ff for A = —1, B = —2, /i = 2, 
z/ = -2. 



order to have i? < with \i < — 4 or ji > — 3, we must have P5 < 0, which means we have 
energy inflow from the bulk. When — 4 < /i < —3, we get P 5 > and energy outflow 
from the brane. Similarly, for A < 0, we get Pb < if z/ > —4. This means that the 
bulk content must have negative pressure. This resembles some sort of five- dimensional 
quintessence field (like a scalar field). For v < —4, the pressure Pb becomes positive 
and we can assume the bulk matter to behave as a regular relativistic fluid. Apparently 
we cannot have both pressures positive, since then we would have two negative powers, 
which is excluded from the constraints obtained above. Notice also that we have managed 
to achieve w e fj — — 1 crossing without the use of any bulk matter which violates energy 
conditions, like phantom fields. The presence of the dark radiation term was also not 
necessary to get the desired phenomenon, which manifests itself as a result of the interplay 
between the bulk pressure and the brane-bulk energy exchange terms. We, thus, have a 
description of dark energy in terms of a bulk energy-momentum tensor, which can drive 
accelerated expansion in the present era and provide a qualitatively correct picture for 
the evolution of the equation of state at late times. 

As a further application of our model, we can show that with proper choice of param- 
eters for the bulk pressure, we can simulate a Chaplygin gas, which is another candidate 
for quintessence. A Chaplygin gas is a perfect four-dimensional fluid which obeys the 
equation of state 

V = -- ■ (28) 
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Assuming energy conservation in the form 



d (pa 3 ) = —pd (a 3 ) , 



(29) 



it is easy to show that for small a, i.e. early times, p oc a~ 3 . This is the behavior we 
expect from dust matter with w = 0. For large a or late times, p = —p, which is the 
equation of state for cosmo logical constant. The Chaplygin gas behaves as dark matter 
at early times and as dark energy (cosmological constant) at late times and is a simple 
example of unified treatment of these two constituents. It is easy to see that the to e // i n 
(123]) can replicate this behavior and thus the bulk content can play a role similar to that 
of a Chaplygin gas. We first want the effective equation of state parameter to be equal 
to —1 at late times (z — 0). This can be satisfied for vA + pB = 0, or 



Additionally, we want w e // — for z » 1. Depending on which is the highest power, 
we must have either p = —3 for v — p > 0, or v = —3 for v — p < 0. The corresponding 
conditions are then: 



For the above choice of parameters, we can reproduce the effects of a Chaplygin gas 
using only bulk matter content. 

Summarizing, we dealt with the problem of simulating the properties of dark energy 
in the context of brane world cosmology. We showed that by assuming the existence of 
an extra dimension and bulk matter, whose energy-momentum tensor is of the form (117j) . 
we can reproduce the observed accelerated expansion of the universe and the equation 
of state for the dark energy. We analyzed in detail the constraints our model parameters 
have to obey in order to get a realistic profile for the evolution of w e ff for small z. The 
crossing of the w e ff = — 1 line occurs naturally, without invoking exotic forms of bulk 
or brane matter that violate energy conditions. Our results also have the advantage of 
not depending on the dark radiation term and they consequently do not fix the sign of 
the C constant. We also showed that this model can reproduce predictions similar to 
those of other dark energy candidates, in particular the Chaplygin gas. We cannot, at the 
moment, fully justify the ansatz (1171) in terms of a phenomenological bulk fluid. As it was 
shown in [13] . this is due to the unknown spatial derivatives of the bulk pressures that 
appear in the equation of energy-momentum conservation for the bulk fluid. A specific 
model for the bulk matter is needed, in order to eliminate this extra freedom. Such a 
realistic model for the bulk matter, which may yield the desired behavior, is currently 
under investigation. 
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-3, z/ = 3f. 
—3, p = 3-g. 
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